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ABSTRACT 

We present a new, robust measurement of the evolving rest-frame ultraviolet (UV) 
galaxy luminosity function (LF) over the key redshift range from z ~ 2 to z ~ 4. 
Our results are based on the high dynamic range provided by combining the Hubble 
Ultra Deep Field (HUDF), CANDELS/GOODS-South, and UltraVISTA/COSMOS 
surveys. We utilise the unparalleled multi-frequency photometry available in this sur¬ 
vey ‘wedding cake’ to compile complete galaxy samples at z ~ 2,3,4 via photometric 
redshifts (calibrated against the latest spectroscopy) rather than colour-colour selec¬ 
tion, and to determine accurate rest-frame UV absolute magnitudes (Afisoo) from 
spectral energy distribution (SED) fitting. Our new determinations of the UV LF ex¬ 
tend from Misoo — —22 (AB mag) down to Mi 5 oo = —14.5, —15.5 and —16 at z ~2, 
3 and 4 respectively (thus reaching ~ 3 — 4 magnitudes fainter than previous blank- 
field studies at z ~ 2 — 3). At z ~ 2 — 3 we find a much shallower faint-end slope 
(a = —1.32 ± 0.03) than reported in some previous studies (a ~ -1.7), and demon¬ 
strate that this new measurement is robust. By z ~ 4 the faint-end slope has steepened 
slightly, to a = —1.43 ± 0.04, and we show that these measurements are consistent 
with the overall evolutionary trend from z = 0 to z = 8. Finally, we find that while 
characteristic number density drops from z ~ 2 to z ~ 4, characteristic luminosity 
(M*) brightens by ~ 1 mag.. This, combined with the new flatter faint-end slopes, has 
the consequence that UV luminosity density (and hence unobscured star-formation 
density) peaks at z ~ 2.5 — 3, when the Universe was ~ 2.5 Gyr old. 

Key words: galaxies: redshifts - galaxies: spectroscopy - galaxies: photometry - 
galaxies: evolution - galaxies: luminosity function 


1 INTRODUCTION 

Ultraviolet (UV) continuum emission provides the most di¬ 
rect tracer of star-formation activity in a galaxy (e.g. Ken- 
nicutt & Evans 2012), albeit it must be corrected for the im¬ 
pact of dust obscuration to derive complete star-formation 
rates. This, coupled with the easy access to the rest-frame 
ultraviolet regime (A = 1500 — 1700 A) provided by opti¬ 
cal observations of galaxies at redshifts z > 1.5, has meant 
that the determination of the evolving UV galaxy luminos¬ 
ity function (LF) has become a key probe of galaxy evolu¬ 
tion and overall cosmic star-formation history (e.g. Bouwens 
et al. 2007; Reddy & Steidel 2009; Robertson et al. 2010; 
McLure et al. 2013; Madau & Dickinson 2014). 

All evidence to date, as provided by a wide range of 
studies discussed later in this paper, indicates that, at red¬ 
shifts 2 ~ 2 — 6, the UV LF is well described by a Schechter 
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function (Schechter 1976). However, as with all LF stud¬ 
ies, the challenge is to assemble galaxy samples of adequate 
size, spanning a large enough dynamic range, and with suf¬ 
ficiently accurate/complete redshift information to robustly 
determine both the bright and faint end of the LF. Thus, 
ideally, large-area surveys are required to adequately sample 
the bright end of the LF (and mitigate the impact of cosmic 
variance) while very deep, small-area surveys are necessary 
to yield the data required to properly constrain the faint-end 
slope, a. This latter quantity is of particular importance if it 
is hoped to perform a reliable luminosity-weighted integral 
of the LF down to faint magnitude limits to determine UV 
luminosity density {puv)- In practice, galaxy selection is also 
a key issue, and it is particularly important to understand 
sample completeness as a function of magnitude and red¬ 
shift, especially if simple colour-colour selection techniques 
are utilized to select samples of UV luminous star-forming 
galaxies. 

After the early pioneering studies of UV luminosity den¬ 
sity indicated it was rising with lookback time out to at least 
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2 ~ 1.5 (e.g. Lilly et al. 1996; Madau et al. 1996; Arnouts 
et al. 2005), the first detailed study of the galaxy UV LF at 
2 ~ 2 — 3 was attempted by Reddy & Steidel (2009). This 
work was based on colour-selection of galaxies from ground- 
based data, and as such was best suited to determining the 
bright-end of the LF at Muv < —18 (AB mag). Nonetheless, 
a key result of this paper was the derivation of an extremely 
steep-faint slope for the UV LF at these intermediate red- 
shifts, with a = —1.73 ± 0.07. This result appears to have 
been confirmed by subsequent, deeper, Hubble Space Tele¬ 
scope (HST) studies following the installation of WFC3 in 
2009; first Oesch et al. (2010) reported a — —1.60 ± 0.21 
at 2 ~ 1.7, and then, very recently, Alavi et al. (2014) used 
WFC3/UVIS combined with the gravitational lensing boost 
provided by the cluster A1689 to probe the UV LF at 2 ~ 2 
down to (lensing corrected) magnitudes Misoo — —13, yield¬ 
ing a = —1.74 ± 0.08. 

However, there has not been universal agreement; for 
example, Hathi et al. (2010) reported a = —1.17 ± 0.40 at 
2 ~ 2.1, and Sawicki (2012) found a = —1.47±0.24 at 2 ~ 2. 
In addition, Weisz, Johnson & Conroy (2014) have recently 
used galactic archaeology techniques to ‘reverse engineer’ 
the stellar populations found in the present-day Local Group 
galaxy population, and conclude in favour of a = —1.35 ± 
0.12 at 2 = 2 - 3. 

Such results are interesting, but this controversy simply 
reinforces the importance of undertaking a new, direct in¬ 
vestigation of the UV LF at this key epoch in cosmic history. 
Moreover, a new measurement of the UV LF at 2 = 2 — 4 is 
timely, given the huge recent improvement in the necessary 
multi-frequency imaging (and supporting spectroscopy) in 
key deep HST and ground-based survey fields. 

In this study we have exploited the combined power of 
the latest optical-infrared data in the Hubble Ultra Deep 
Field (HUDF) (reaching ~ 29.5 mag over 4.5arcmin^), 
the CANDELS/GOODS-S field (reaching ~ 27.5 mag over 
170arcmin^), and the UltraVISTA/COSMOS field (reach¬ 
ing ~ 26mag over ~ 0.5deg^). In recent years, driven by 
the rapid improvements in the near-infrared depth arising 
from the advent of WFG3/IR and VISTA imaging (e.g. Gro- 
gin et al. 2011; McGracken et al. 2012; Koekemoer et al. 
2013) the unparalleled data in these key fields have been 
very actively exploited in the study of the UV LF at higher 
redshifts, 2 = 5 — 8 (e.g. McLure et al. 2010, 2013; Oesch 
et al. 2010, 2013; Bouwens et al. 2011, 2015; Finkelstein 
et al. 2010, 2015; Dunlop 2013; Bowler et al. 2014, 2015). 
However, the full multi-frequency datasets have not recently 
been properly applied to revisit the measurement of the UV 
LF in the redshift range 2 = 2 — 4. Onr new study thus aims 
to rectify this situation, and to also take advantage of new 
optical and near-infrared spectroscopy (including WFG3/IR 
grism spectroscopy; Skelton et al. 2014; Morris et al. 2015) to 
help produce the most reliable photometric redshifts (crucial 
both for robust galaxy sample selection, and for the accurate 
determination of Misoo or M 1700 for each source). 

The remainder of this paper is structured as follows. In 
Section 2 we summarize the available imaging/photometric 
data for our 3-tier ‘wedding-cake’ survey. Next, in Section 3, 
we summarize the available spectrosocopic data, and explain 
how we have used robust photometric redshift estimation to 
assemble a combined sample of ~ 35, 000 galaxies in the red¬ 
shift range of interest, 1.5 < 2 < 4.5. We then proceed, in 


Section 4, to analyse the resulting galaxy dataset to derive 
the UV LF, focussing first on the detailed shape (especially 
the faint-end slope, a) of the LF at 2 ~ 2, and then expand¬ 
ing the analysis up to 2 ~ 4. In Section 5 we compare our 
results to those of other authors, place our findings in the 
context of published results on the evolution of the UV LF 
up to 2 ~ 8, and discuss the implications of our derived UV 
luminosity density, puv, over the redshift range 2 = 2 — 4. 
Our conclusions are summarized in Section 6. 

Throughout we assume a flat cosmology with Do = 0.3, 
Da = 0.7 and Ho = 70kms“^Mpc“^, and give all magni¬ 
tudes in the AB system (Oke 1974; Oke & Gunn 1983). 


2 IMAGING DATA AND PHOTOMETRY 

In this section we summarize the properties of the three 
surveys used in this study, and explain how the galaxies were 
selected, and their multi-frequency photometry measured. 

Throughout we refer to the HST AGS and WFG3/IR 
filters F435W, F606W, F775W, F814W, F850LP, F098M, 
F105W, F125W, F140W and F160W as B 435 , Ueoe, * 775 , 
*814, 2850, Yo98, Ylos, Ji25, J^i40 and Hloo respectively, the 
VLT VIMOS-U and Hawk-I/K-215 filters as U and Ks re¬ 
spectively, the GFHT MegaGam optical filters as u, g, r, 
i', z\ the reddest Subaru Suprime-Gam filter as 2 , the four 
broad-band VISTA near-infrared filters as Y, J, H and Ks, 
and the Spitzer IRAG first two channels as IRACo.otim and 

IRACi.Ofim- 

The photometric depths of the imaging data in each 
filter, for each field, are summarized for convenience in Table 
1, with further details for each field given in the following 
subsections. 

2.1 HUDF 

The deepest dataset we analyse, which is crucial for con¬ 
straining the faint end of the galaxy UV LF, is the multi¬ 
band imaging of the Hubble Ultra Deep Field completed in 
2012 (HUDF12), covering an area of ~ 4.5arcmin^. This 
latest comprehensive dataset consists of the deepest near- 
infrared imaging obtained with HST WFG3/IR from the 
HUDF09 and HUDF12 programmes (Bouwens et al. 2010; 
Ellis et al. 2013; Koekemoer et al. 2013; Dunlop et al. 2013; 
Robertson et al. 2013), and the original optical HST AGS 
imaging (Beckwith et al. 2006) supplemented by new deep 
*814 imaging. To maximise wavelength coverage, we have 
supplemented the HST data with the public VLT GOODS- 
S U imaging (Nonino et al. 2009), the Hawk-I Ks imaging 
from the HUGS programme (Fontana et al. 2014), and the 
deepest available IRAG imaging (McLure et al. 2011; Ashby 
et al. 2013). 

Galaxy detection and photometry from this deep HST 
imaging dataset was undertaken using SEXTRACTOR v2.8.6 
(Bertin & Arnouts 1996) in dual-image mode with Hioo 
as the detection image and the flux-iso as the observed 
isophotal flux. In order to obtain consistent resolution- 
matched photometry, the lower resolution U, Ks and IRAG 
images were deconfused using the technique described in 
McLure et al. (2011). 

Our final photometric catalogue for the HUDF includes 
the photometry of 2864 sources with iJieo < 29.5, measured 
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Table 1. A summary of the photometry used in this study, giving 
the 5 -(t detection limits in each filter/field as appropriate. The 
filter names are as summarized at the beginning of Section 2. For 
the HST photometry the depths given refer to total magnitudes, 
as derived from small-aperture magnitudes assuming point-source 
corrections (see McLure et al. 2013). The Spitzer IRAC depths 
also refer to total magnitudes, as derived using TPHOT (Merlin 
et al. 2015). The depths for the ground-based photometry are 
based on 2-arcsec diameter aperture measurements. 


Filter 


Survey 



HUDF 

CANDELS/ 

GOODS-S 

UltraVISTA/ 

COSMOS 

U 

28.0 

28.0 


U 



27.0 

B 435 

29.7 

28.0 


9 



27.1 

Veoe 

30.2 

28.4 


r 



26.6 

*775 

29.9 

27.8 


i' 



26.3 

^814 

29.8 



z' 



25.4 

z 



26.4 

^850 

29.1 

27.5 


Ul05 

29.7 

27.9 


Y 



25.1 

J 125 

29.2 

27.7 


J 



24.9 

Jl40 

29.2 



H 16 O 

29.2 

27.3 


H 



24.6 

Ks 

26.5 

26.5 

24.8 

IRAC 3.6 

26.5 

26.5 

25.2 

IRAC 4.8 

26.3 

26.3 

25.2 


in the U, B 435 , Veoe, irrs, *8i4, zsso, ^loS) >^ 125 , Jno, Hiqq, 
and IR^C 4 :. 5 ^m bands. 

2.2 CANDELS/GOODS-S 

To provide the next tier of the survey ‘wedding cake’, we 
have used the publicly-available HST WFC3/IR and HST 
ACS imaging of the Great Observatories Origins Deep Sur¬ 
vey South (GOODS-S) held provided by the Cosmic Assem¬ 
bly Near-Infrared Deep Extragalactic Legacy Survey (CAN- 
DELS) (Grogin et al. 2011; Koekemoer et al. 2011; Wind¬ 
horst et al. 2011), and the associated pre-existing HST opti¬ 
cal (Giavalisco et al. 2004; Riess et al. 2007), ground-based 
VLT IZ-band (Nonino et al. 2009) and As-band (Retzlaff et 
al. 2010; Fontana et al. 2014), and Spitzer IRAC imaging 
(Ashby et al. 2013), as summarized by Guo et al. (2013). 

Consistent with the production of the HUDF catalogue, 
the sources were detected, and their isophotal huxes in the 
HST bands measured using SEXTRACTOR v2.8.6 (Bertin & 
Arnouts 1996) in dual-image mode, again with iLieo as the 
detection image. As described in Guo et al. (2013), in this 
held the Template FITting (tfit) method (Laidler et al. 
2007) has been applied to generate the matched photometry 
from the lower angular resolution U, Ks and IRAC imag¬ 
ing. The GOODS-S catalogue provided by Guo et al. (2013) 
contains 34930 sources in an area of 173arcmin^, with pho¬ 


tometry in the U, 5435, Veoe, * 775 , .^850, logs, Lios, J 125 , 
H 160 , Ks, IRACiSiim and IRACi.e,^^ bands. 

2.3 UltraVISTA/COSMOS 

Data Release 2 (DR2) of the Ultra VISTA survey (Mc¬ 
Cracken et al. 2012) provides deep near-infrared imaging 
in 4 deep strips which overlap ~ 0.7 deg^ of the area also 
covered by the HST imaging of the Cosmological Evolution 
Survey (COSMOS; Scoville et al. 2007), and by the Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS) Mega- 
Cam deep optical imaging. As discussed in Bowler et al. 
(2012, 2014, 2015), for the central square degree covered by 
the CFHTLS D2 imaging, very deep Subaru Suprime-Cam 
2 -band imaging has also been obtained (Furusawa et al., in 
preparation). 

The photometry of the sources in the UltraV- 
ISTA-1-COSMOS imaging has been measured in 2" diameter 
apertures using SEXTRACTOR in dual image mode. The de¬ 
tection image used in this case is the i'—band image from the 
T0007 release of the CFHTLS. These fluxes have been con¬ 
verted to total using the i'—band FLUX_AUTO parameter, 
and in this case the new tphot code (Merlin et al. 2015) has 
been used to obtain deconfused, resolution-matched IRAC 
fluxes. 

The final UltraVISTA-|-COSMOS catalogue utilised 
here contains 89614 galaxies with i! < 26, selected from 
an area of 0.482 deg^ (reduced to an effective survey area of 
0.292 deg^ after masking for bright objects, and diffraction 
spikes etc), and provides photometry in the u, g, r, i', and 
2 '—bands from the CFHTLS, the 2 —band from Suprime- 
Cam on Subaru, the Y, J, H, and Ks bands from Ultra- 
VISTA, and deconfused IRACs.e/im and IRAC 4 . 5 ,im pho¬ 
tometry from a combination of the Spitzer Extended Deep 
Survey (SEDS; PI: Fazio; Ashby et al. 2013) and the Spitzer 
Large Area Survey with Hyper-SuprimeCam (SPLASH; PI: 
Capak). 


3 REDSHIFT INFORMATION 

In an ideal world, the galaxy LFs would be derived from 
large, complete, galaxy samples with complete spectroscopic 
redshift information. However, at the depths of interest 
here, even semi-complete spectroscopic redshift information 
is clearly (currently) not possible. Consequently, the robust¬ 
ness of our final results depends crucially on the reliability 
and accuracy of photometric redshifts derived for the galax¬ 
ies uncovered for each of the three survey fields discussed 
above. 

To test/optimize the derived photometric redshifts, sub¬ 
sets of objects with reliable spectroscopic redshifts are re¬ 
quired. We have therefore assembled the latest catalogues 
of robust spectroscopic redshifts in the HUDF, GOODS-S, 
and COSMOS fields. The results of this search are summa¬ 
rized in the hrst sub-section below. 

We then explain the steps taken to optimize photomet¬ 
ric redshifts, the final procedure adopted, and quantify the 
reliability and accuracy of our final photometric redshifts via 
comparison with the spectroscopic database (with further 
details and comparisons with other published photometric 
redshift catalogues given in Appendix A). 
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Figure 1. The distribution of the 2799 high-quality spectroscopic 
redshifts in the latest spectroscopic sample we have assembled for 
the CANDELS/GOODS-S field. The red histogram indicates the 
new WFC3/IR near-infrared grism redshifts as determined by the 
3D-HST survey (Skelton et al. 2014) and CANDELS (Morris et 
al. 2015), which have helped significantly to fill in the ‘redshift 
desert’. Of the 2799 objects with high-quality spectroscopic red¬ 
shifts shown here, 218 lie within the HUDF. 

In the third subsection below we describe our final 
combined sample of galaxies with redshifts in the range 
1.5 < « < 4.5. This sample is used in the LF analyses pre¬ 
sented in the remainder of this paper. 


3.1 Spectroscopic Redshifts 

Despite the fact that the HUDF, GOODS-S and COSMOS 
fields have been targeted by several of the most dense and 
deep spectrosopic surveys ever undertaken, typically at most 
~ 10% of the galaxies in our photometric samples possess 
high-quality spectroscopic redshifts. Nonetheless, spectro¬ 
scopic redshift information is crucial for refining and quan¬ 
tifying the accuracy (e.g. a) and reliability (e.g. number of 
catastrophic outliers) of photometric redshifts. 

We have therefore assembled catalogues including the 
very latest spectroscopic redshift information in each field. 
Within the GOODS-S field (including the HUDF) we have 
assembled a sample of 2799 galaxies with high-quality red¬ 
shifts (218 of which lie within the area covered by the 
WFC3/IR imaging of the HUDF). We have confined our se¬ 
lection to only the very highest quality flags assigned to the 
redshifts obtained by each study in the literature, and the re¬ 
sulting redshift distribution of our final spectroscopic galaxy 
sub-sample (after removal of any stars or AGN) is shown in 
Fig. 1. As illustrated by the blue histogram in Fig. 1, the ma¬ 
jority of this spectroscopic redshift information (1917 red¬ 
shifts) has been obtained from ground-based optical spec¬ 
troscopy (Balestra et al. 2010; Cimatti et al. 2008; Gristiani 
et al. 2000; Groom et al. 2001; Dohetry et al. 2005; Szokoly 


et al. 2004; Le Fevre et al. 2004; Mignoli et al. 2005; Roche et 
al. 2006; Silverman et al. 2010; Strolger et al. 2004; Vanzella 
et al. 2008). However, as illustrated by the red histogram 
in Fig. 1, it can also been seen that recent HST WFG3/IR 
near-infrared grism spectroscopy (Skelton et al. 2014; Morris 
et al. 2015) has now made an important contribution to the 
redshift coverage in this deep field (982 redshifts), in partic¬ 
ular helping to fill in the traditional ‘redshift desert’ between 
2 ~ 1.2 and 2 ~ 2, where relatively few strong emission lines 
are accessible in the optical regime. 

Within the COSMOS field we utilized a sample of 1877 
high-quality redshifts as provided by the public 2 -GOSMOS 
survey (Lilly et al. 2007). 

3.2 Photometric Redshifts 

3.2.1 Method 

To determine photometric redshifts we used the pub¬ 
lic galaxy template-fitting code Le Phare[^ (PHotometric 
Analysis for Redshift Estimate; Ilbert et al. 2006). To en¬ 
sure the proper treatment of weak/non-detections, we fit¬ 
ted in flux-density rather than magnitude space. To ac¬ 
count for dust obscuration/reddening we assumed the dust- 
attenuation law of Calzetti et al. (2000), allowing reddening 
to vary over the range 0.0 < E{B — V) < 0.5 in steps of 
AE{B — V) = 0.1. We also included IGM absorption as¬ 
suming the models of Madau (1995). 

For each field we proceeded in four stages. First, to 
avoid too much weight being placed on individual photo¬ 
metric detections, and to allow for remaining low-level sys¬ 
tematic errors, we set a minimum of error of 3% on all optical 
and near-infrared photometry, and a minimum error of 10% 
on all IRAC photometry. Next, we utilized the galaxy SED 
templates provided by the evolutionary synthesis models of 
Bruzual & Chariot (2003) (BC03), without emission lines, 
and adjusted the photometric zero-points until the accu¬ 
racy of the photometric redshifts was maximised (as judged 
by comparison with the high-quality spectroscopic redshifts 
discussed above). After this we explored the use of a range 
of different galaxy templates, before determining that the 
the PEGASEv2.0 models (Fioc & Rocca-Volmerange 1999), 
with emission lines switched on, produced the most accu¬ 
rate final photometric redshifts. Finally, we determined the 
photometric redshift for each galaxy by searching the red¬ 
shift range 2 = 0 — 10, and distinguished between acceptable 
and unacceptable photometric redshifts based on an analysis 
of the distribution of minimum (resulting in the accep¬ 
tance of final SED fits with minimum < 50 in the fields 
with HST photometry, and minimum x^ < 20 in the Ultra- 
VISTA/COSMOS field). This level of quality control led to 
the exclusion of < 5% of galaxies in the photometric samples 
from the final sample with trusted photometric redshifts. 

Below we summarize the results obtained in each field. 


^ http://www.cfht.hawaii.edu/-arnouts/LEPHARE/lephare. 
html 
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Figure 2. The comparison of our new photometric redshifts ver¬ 
sus the high-quality spectroscopic redshifts in each of the three 
survey fields. The top plot for the HUDF contains 210 galaxies, 
the central plot for CANDELS/GOODS-S contains 2677 galaxies, 
while the bottom plot for UltraVISTA/COSMOS contains 1671 
galaxies. As summarized by the statistics in each panel, the ac¬ 
curacy of our photometric redshifts is comparable with the very 
best ever achieved for high-redshift galaxy surveys; see Section 3 
and Appendix A for futher details. 





Figure 3. A second comparison of our new photometric redshifts 
versus the high-quality spectroscopic redshifts in each of the three 
survey fields, this time confined to 2 : > 1.5. The HUDF plot con¬ 
tains 72 galaxies, the CANDELS/GOODS-S plot contains 668 
galaxies, while the UltraVISTA/COSMOS plot contains 71 galax¬ 
ies. As summarized by the statistics in each panel, the accuracy 
of our photometric redshifts remains comparable with the very 
best ever achieved for high-redshift galaxy surveys; see Section 3 
and Appendix A for futher details. 
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3.2.2 HUDF photometric redshifts 

As outlined above, we refined the HUDF photometric zero- 
points by fitting the photometry with the BC03 models; the 
derived zero-point offsets were all smaller than 0.1 mag. Af¬ 
ter application of Le Phare with emission lines, we derived 
acceptable (x^ < 50) photometric redshifts for 2730 galaxies 
with H < 29.5 in the HUDF. 

The accuracy of the derived photometric redshifts in 
the HUDF is illustrated in the top panel of Fig. 2, which 
shows Zphot V Zapec for the 210 sources in the field with 
secure spectroscopic redshifts, and acceptable photometric 
redshifts. The outlier fraction is 4.2%, and (tnmad = 0.026. 
For redshifts confined to 2 > 1.5 the Zphot v Zapac plot is 
shown in the top panel of Fig. 3 (see Appendix A for further 
details). 

3.2.3 CANDELS/GOODS-S photometric redshifts 

To maximise the reliability of the photometric redshifts 
for this sample, we confined our attention to sources with 
Hieo < 26 in the CANDELS Wide region, and to sources 
with Hiqo < 27 in the CANDELS Deep region (which covers 
the central ~ 55arcmin^ of GOODS-S), and then again re¬ 
fined the photometric zeropoints by fitting the photometry 
with the BC03 models; the derived zero-point offsets were 
again all smaller than 0.1 mag. 

After application of Le Phare with emission lines, we 
derived acceptable {p < 50) photometric redshifts for 10987 
galaxies with iL < 27 in the GOODS-S Deep field, and for 
27460 galaxies with iL < 26 in the GOODS-S Wide field. 

The accuracy of the derived photometric redshifts in 
CANDELS/GOODS-S is illustrated in the middle panel of 
Fig. 2, which shows Zphot v Zapac for the 2677 sources in 
the field with secure spectroscopic redshifts, and accept¬ 
able photometric redshifts. The outlier fraction is 3%, and 
(Tnmad = 0.027. For redshifts confined to z > 1.5 the Zphot 

V Zapac plot is shown in the middle panel of Fig. 3 (see Ap¬ 
pendix A for further details). 

3 . 2.4 UltraVISTA/COSMOS photometric redshifts 

As with the HST-hased catalogues, initially the photomet¬ 
ric redshifts were computed by fitting the ground-based and 
Spitzer photometry to the BG03 models, in order to adjust 
the photometric zero-points until maximum redshift accu¬ 
racy was achieved; all the derived zero-point offsets were 
again all smaller than 0.1 mag. 

For consistency with the HST HUDF-I-GOODS-S data 
analysis, the photometric redshifts and rest-frame abso¬ 
lute UV magnitudes were then again recomputed using Le 
Phare with emission lines, yielding acceptable (x^ < 20) 
photometric redshifts for 88789 galaxies with i' < 26 in the 
UltraVISTA/GOSMOS field. 

The accuracy of the derived photometric redshifts in 
Ultra VISTA/GOSMOS is illustrated in the bottom panel of 
Fig. 2, which shows Zphot v Zapac for the 1671 galaxies in 
the field with secure spectroscopic redshifts, and acceptable 
photometric redshifts. The outlier fraction is only 2.1%, and 
(Tnmad = 0.026. For redshifts confined to 2 > 1.5 the Zphot 

V Zapaa plot is shown in the bottom panel of Fig. 3 (see Ap¬ 
pendix A for further details). 



Figure 4. The distribution of photometric redshifts for our final 
galaxy sample comprising a total of 36051 galaxies in the redshift 
range 1.5 < Zphot < 4.5. The different coloured histograms show 
the redshift distribution subdivided by survey-field/depth, with 
1549 galaxies from the HUDF with ifieo < 29.5 (red), 4465 galax¬ 
ies from CANDELS/GOODS-S Deep with H\qq < 27 (cyan), 
6947 galaxies from CANDELS/GOODS-S Wide with H\eo < 26 
(blue), and 23090 galaxies from the Ultra VISTA/COSMOS field 
with i' < 26.0 (yellow). 


3.3 Final Galaxy Sample 

The final galaxy sample consists of 36051 galaxies selected 
to he in the redshift range 1.5 < 2 < 4.5, and consists 
of: i) 1549 galaxies from the HUDF with Hieo < 29.5, 
ii) 11412 galaxies from GOODS-S (comprising 4465 from 
CANDELS/GOODS-S Deep with Hieo < 27 and 6947 from 
CANDELS/GOODS-S Wide with Hmo < 26), and iii) 
23090 galaxies from the UltraVISTA/GOSMOS field with 
i' < 26.0. 

The redshift distributions of these final HUDE, 
GANDELS/GOODS-S and UltraVISTA/GOSMOS samples 
are shown in Fig. 4. 

Absolute magnitudes at UV rest-frame wavelengths 
\raat = 1500 A and \raat = 1700 A for use in the subse¬ 
quent LF analyses were computed from the Le Phare SED 
fits using a 100 A top-hat synthetic filter centred at the ap¬ 
propriate rest wavelength. 


4 THE GALAXY LUMINOSITY FUNCTION 

Armed with redshifts and absolute UV magnitudes for over 
35,000 galaxies in the redshift range 1.5 < z < 4.5, we can 
now derive the rest-frame UV galaxy LF, exploring its form 
and evolution from 2 = 2 to 2 = 4. To aid comparison 
with previous results in the literature, we derive LFs at both 
Xraat = 1500 A and Xreat = 1700 A as required, but for our 
final calculations of the evolution of LF parameters, and the 
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evolution of luminosity density (and hence star-formation 
rate density), we focus on 1500 A at all redshifts. 

In the first subsection below we outline the (straight¬ 
forward) method we have adopted to determine the non- 
parametric binned form of the LF, and then the parametric 
form (in this case the Schechter function). 

We then present and discuss our results at various red- 
shifts, in part to facilitate comparison with the literature. 
Specifically, we consider first the UV LF at 1500 A in the 
redshift range 1.5 < a < 2.5 {z ~ 2) (in particular focussing 
on the faint-end slope, a), before considering separately the 
1500 A LF at 1.5 < z < 2 ( 2 ; ~ 1.7) and the 1700 A LF at 
2 < z < 2.5 {z ~ 2.2). We then move on to determine the 
evolution of the 1500 A galaxy LF over the redshift range 
2 ~ 1.5 — 4.5, in three redshift bins of width Az = 1 (i.e. 
corresponding to 2 ~ 2, z ~ 3 and 2 ~ 4). 


4.1 Method 

Various techniques can be ultilized to derive the LF, but 
here we have sufficiently extensive and dense coverage of 
the luminosity-redshift plane to obtain a non-parametric 
estimate by applying the straightforward Vmax estimator 
(Schmidt 1968) given by: 

<).(M)dM = ^[—-] (1) 

i 

where the sum is over all galaxies in the given redshift and 
absolute magnitude bin (chosen here to have a width of 
AM = 0.5mag), Vmax for each galaxy is set by the up¬ 
per redshift limit of the bin unless the source drops out of 
the sample before that redshift is achieved, and C is the 
completeness factor for each source. The errors on the de¬ 
rived number density in each bin are here assumed to be 
Poissonian. 

The completeness factor corrects for incompleteness 
caused both by the fact that significant regions of the imag¬ 
ing are in practice inaccessible for high-redshift object se¬ 
lection (i.e. areas masked due to the presence of bright fore¬ 
ground galaxies or stars/diffraction spikes) and by photo¬ 
metric scatter (which obviously impacts most seriously on 
the faintest magnitude bins). This has been calculated by 
Monte Carlo source injection and retrieval simulations, and 
over most of the magnitude range in each sample transpires 
to be ~ 70% in the HUDF and CANDELS/GOODS-S fields, 
and ~ 60% in the UltraVISTA/COSMOS field. 

Obviously, within each of the three survey fields uti¬ 
lized here, incompleteness becomes more serious as the de¬ 
tection limit is approached, and the impact of photometric 
scatter becomes significant. However, in the present study 
the impact of this is minimal, as there is sufficient overlap 
between the regions of the luminosity-redshift plane cov¬ 
ered by the different surveys that we can, for example, dis¬ 
card all seriously incomplete faint bins from the UltraV¬ 
ISTA/COSMOS survey in favour of the first well-sampled 
brighter bins from the CANDELS/GOODS-S survey (and 
similarly ensuring the LF determination is dictated by the 
HUDF before CANDELS/GOODS-S becomes seriously in¬ 
complete). At the very faint end we neglect all bins delivered 
by the HUDF in which photometric scatter results in a com¬ 
pleteness < 90%. 




Figure 5. The galaxy rest-frame UV LF at 2 ~ 2. The upper 
panel shows the new 1500 A LF as derived from our combined 
galaxy sample in the redshift range 1.5 < 2 < 2.5. The blue circles 
with errors indicate the number densities from the Vmax estima¬ 
tor (see Table 2) and the blue line is our best-fitting Schechter 
function. The red and pink lines are the Schechter functions at 
2 ~ 2 reported by Oesch et al. (2010) and Alavi et al. (2014) re¬ 
spectively, both of whom derived a much steeper faint-end slope 
at 2 ~ 2. The orange line shows the 2 ~ 2 LF as inferred by 
Weisz et al. (2014); this has a significantly shallower faint-end 
slope, in excellent agreement with the value of a deduced from 
our new determination (see Section 4.2.1). The lower panel shows 
how the fitted value of a depends on the limiting absolute magni¬ 
tude down to which the fitting is performed. It can be seen that 
the derived faint-end slope stabilises at M1500 > —17, settling to 
a secure and robust value of q = —1.32 ± 0.03. 


Finally, having derived the non-parametric LF from the 
combination of our three surveys, we fit the binned values 
with a Schechter function (Schechter 1976): 

</(M) = (2) 


© 2014 RAS, MNRAS 000, [^?? 





8 Parsa et al. 


where (j)*, M* and a indicate respectively the normalisa¬ 
tion coefficient, the characteristic magnitude and the faint 
end slope of the LF, and derive confidence intervals on the 
parameters. 

4.2 The galaxy UV LF at z ~ 2 

4.2.1 1.5 <z <2.5 

We first derive a new measurement of the galaxy rest-frame 
UV LF at 2 ~ 2, based on all galaxies in our combined sam¬ 
ple with photometric redshifts in the range 1.5 < z < 2.5. 
Fig. 5 shows the resulting LF at Arest = 1500 A, including 
our best-fitting Schechter function. Here, bins brighter than 
Misqo = —20 are dominated by the UltraVISTA COSMOS 
sample, while at the faintest magnitudes the HUDF sample 
allows us to extend the UV LF down to M 1500 — —14.5, 
which is ~ 5 mag fainter than achieved by Oesch et al. 
(2010). The extension of the 2 : ~ 2 LF to comparably faint 
magnitudes has only previously been reported by Alavi et 
al. (2014). However, this was only achieved with the aid of 
the gravitational leasing provided by the cluster Abell 1689, 
resulting in very small effective survey volumes and conse¬ 
quently much poorer S/N than achieved here. 

The most striking result of our new 2 ~ 2 LF deter¬ 
mination, as shown in the upper panel of Fig. 5, is that we 
hnd a much shallower faint-end slope (a = —1.32 ± 0.03) 
than reported by either Oesch et al. (2010) or Alavi et al. 
(2014) who found a = —1.60 ± 0.51 and a = —1.74 ± 0.08 
respectively. The lower panel in Fig. 5 shows that an accu¬ 
rate measurement of a requires good sampling of the LF at 
magnitudes fainter than M 1500 — —17, after which the fitted 
value of a stabilises and yields a robust measurement. It is 
therefore not surprising that Oesch et al. (2010) deduced an 
erroneously steep faint-end slope given the limited depth of 
the data utilised in that study. 

Our best-fitting values of the other Schechter parame¬ 
ters at 2 ~ 2 are M* = —19.68 ± 0.05 and p = 7.02 ± 0.66 
(x 10“®Mpc“®mag“'^). Also shown in Fig. 5 is the 2 ~ 2 UV 
LF inferred by Weisz et al. (2014) from local-group galac¬ 
tic archaeology. Interestingly it is this ‘reverse engineered’ 
LF which agrees best with our new direct determination, as 
Weisz et al. (2014) also infer a very similar, shallow faint-end 
slope of a = —1.36 ± 0.11. 

4 . 2.2 1.5 < 2 < 2 and 2 < 2 < 2.5 

To further facilitate comparison with previous studies, we 
next compute the UV LF within finer redshift bins, 1.5 < 
2 < 2 and 2 < 2 < 2.5. We calculate the LF in the latter 
redshift bin at Xrest = 1700 A to simplify direct compari¬ 
son with the results of Reddy & Steidel (2009) and Sawicki 
(2012), both of whom calculated the LF at 1700 A based 
on colour selection sampling an effective redshift window 
1.9 < 2 < 2.7. 

Our results at 2 ~ 1.7 and 2 ~ 2.2 are shown in Fig. 6. 
As in Fig. 5 the upper panels show the binned LF and best- 
htting Schechter function, while the lower panels show the 
derived value of faint-end slope, a, as a function of the ab¬ 
solute magnitude limit down to which the fitting is per¬ 
formed. Again it can be seen that the derived value of a 
only stabilises at Muv > —17, and that studies reaching 


only Muv — —19 are likely to yield an eroneously steep 
faint-end slope. 

Our derived Schechter parameter values for the 1500 A 
LF at 2 ~ 1.7 are M* = -19.61 ± 0.07, p = 6.81 ± 
0.81 (x 10“®Mpc“®mag“^), and a = —1.33 ± 0.03. In the 
left-hand panel of Fig. 6 we also show the Schechter function 
derived by Oesch et al. (2010) in the same redshift inter¬ 
val; it can be seen that while the bright end is comparable, 
Oesch et al. (2010) inferred a much steeper faint-end slope 
of a = -1.6 ±0.21. 

For the 1700 A LF at 2 ~ 2.2 we find M* = —19.99 ± 
0.08, p = 6.20 ± 0.77(xl0"^Mpc"^mag"^), and a = 

— 1.31 ± 0.04. In the right-hand panel of Fig. 6 we also 
show the Schechter functions derived by Reddy & Steidel 
(2009) and Sawicki (2012) from very similar redshift ranges. 
While our LF matches that derived by Reddy & Steidel 
(2009) around the break, the faint-end slope derived by 
Reddy & Steidel (2009) was clearly much steeper, with a = 

— 1.73 ±0.07. The faint-end slope derived by Sawicki (2012) 
was somewhat shallower (although still steeper than our new 
derivation), but rather uncertain {a = —1.47 ± 0.24). More¬ 
over, it can also be seen that the bright-end of the LF as 
derived by Sawicki (2012) also deviates significantly from 
our new results. 


4.3 The galaxy UV LF at z ~ 3 and z ~ 4 

We now extend our study of the galaxy UV LF out to higher 
redshift, considering also the redshift bins 2.5 < 2 < 3.5 and 
3.5 < 2 < 4.5, in order to explore how the LF evolves over 
the crucial redshift range 2 ~ 2 — 4. 

We have focussed on 1500 A and our results are shown 
in Fig. 7, alongside the 2 ~ 2 1500 A LF which was shown in 
Fig. 5. The binned number densities derived from the Vmax 
method at 2 ~ 2, 3 and 4, as shown in Fig. 7, are tabulated in 
Table 1, and our derived best-fitting Schechter function pa¬ 
rameter values at all three redshifts are included in Table 2 
(along with various values from the literature, as discussed 
further in Section 5). 

Again, for comparison, in Fig. 7 we also overplot other 
recent determinations of the 1500 A LF at these redshifts. 
Our derived LFs at 2 ~ 3 and 2 ~ 4 appear to agree reason¬ 
ably well with previous measurements around the break lu¬ 
minosity but, as at 2 ~ 2, we find a shallower faint-end slope, 
a, than most previous studies; again, we agree best with the 
results inferred from the local Universe by Weisz et al. (2014) 
(although we note that, at 2 ~ 3, Weisz et al. (2014) have 
adopted the bright end of the 1500 A LF given by Reddy 
& Steidel (2009)). At 2 ~ 4 we hnd a shallower faint-end 
slope than most previous studies (i.e. a = —1.43 ± 0.04), 
but (j>* and M* are in excellent agreement with the results 
of Bouwens et al. (2014) (see Table 3). 

In Fig. 8 we plot our derived Schechter parameter val¬ 
ues, with 1-CT and 2-g single-parameter conhdence intervals, 
for 2 ~ 2, 2 ~ 3 and 2 ~ 4, while in Fig. 9 we overplot 
the 1500 A LFs at these three redshifts. These plots high¬ 
light evolutionary trends in the UV LF over this key redshift 
range, which we discuss further below in Section 5. 
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Figure 6. Upper-Left: The 1500 A LF as derived for galaxies in the redshift bin 1.5 < z < 2. The blue data points show the number 
densities measured using the Vmax method and the solid blue line is our best-fitting Schechter function. The yellow solid line with 
the steeper faint-end slope is the best-fitting Schechter LF measured by Oesch et al. (2010) over the same photometric redshift range. 
Upper-Right: The 1700 A LF derived over the photometric redshift range 2 < z < 2.5. Again the blue data points show the binned 
number densities as derived from the Vmax method, while the solid blue line is our best-fitting Schechter function. The pink and yellow 
lines indicate, respectively, the Schechter function LFs derived by Reddy & Steidel (2009) and Sawicki (2012), based on colour-selection. 
The lower panels again show how the fitted value of a depends on the limiting absolute magnitude down to which the fitting is performed. 


5 DISCUSSION AND IMPLICATIONS 


5.1 The evolution of the LF from z ~ 2 to z ~ 4 

As can be seen in Fig. 8 and Fig. 9, while the LF displays 
relatively little evolution between 2 ; ~ 2 and 2 ~ 3, there is a 
clear drop in (j)* (by a factor ~ 2.5) between 2 ~ 3 and 2 ~ 4. 
We also find, however, that M* brightens steadily over this 
redshift range, by ~ 1 mag.. In terms of luminosity density, 
this brightening more than offsets the decline in cj>* up until 
2 ~ 3, but by 2 ~ 4 the more dramatic drop in 0* dominates 


the evolution, and luminosity density undoubtedly declines 
significantly by 2 ~ 4. 

The modest evolution in cj)* and the ~ 0.5 mag. bright¬ 
ening in M* seen between 2 ~ 2 and 2 ~ 3 is very similar to 
the evolution reported by Reddy & Steidel (2009), with the 
main difference being that our (j>* values are systematically 
higher, and our M* values systematically dimmer due (at 
least in part) to our significantly shallower best-fitting faint- 
end slope a. This is made clear in Fig. 10, where we place our 
results in the context of several other recent studies. Here 
it can be seen that our inferred values of <))* at 2 ~ 2 and 
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Figure 7. Our new rest-frame UV (1500 A) galaxy luminosity functions at z ~ 2, 3 and 4. The (blue, green, red) data points indicate 
the values derived via the Vmax estimator, with the colour-matched solid lines showing the best-fitting Schechter functions. The values 
corresponding to the data points and their errors are tabulated in Table 2, while the values of the best-fitting Schechter parameters are 
given in Table 3, and plotted in Fig. 8. The vertical dashed line in each panel at M 1500 = is shown simply to indicate the typical 
absolute magnitude limit reached by previous studies at these redshifts, while the background shading indicates the absolute magnitude 
regimes in which the three different surveys (UltraVISTA/COSMOS, CANDELS/GOODS-S, HUDF) make the dominant contribution 
to our new measurement of the LF at each redshift. For comparison purposes, at each redshift we also show several luminosity functions 
from the literature as indicated in the legend, and discussed in Section 4.3. 


Table 2. The rest-frame UV (1500 A) galaxy luminosity functions at 2 : ~ 2, 3 and 4, as measured via the Vmax estimator; these values 
are plotted in the three panels of Fig. 7, and in Fig. 9. 


Misoo ~ 2) 0 ( 2 : ~ 3) 

/Mpc“^mag~^ /Mpc~^mag“^ 


0 (z ~ 4 ) 

/Mpc“®mag“^ 


-23 

- 



- 



0.000001 

± 

0.000000 

-22.5 

- 



0.000003 

± 

0.000001 

0.000010 

± 

0.000002 

-22 

0.000012 

± 

0.000027 

0.000023 

± 

0.000004 

0.000043 

± 

0.000005 

-21.5 

0.000034 

± 

0.000045 

0.000117 

± 

0.000008 

0.000154 

± 

0.000010 

-21 

0.000152 

± 

0.000094 

0.000462 

± 

0.000016 

0.000475 

± 

0.000017 

-20.5 

0.000555 

± 

0.000181 

0.001462 

± 

0.000107 

0.001087 

± 

0.000096 

-20 

0.001654 

± 

0.000124 

0.002511 

± 

0.000140 

0.001709 

± 

0.000120 

-19.5 

0.003467 

± 

0.000165 

0.003830 

± 

0.000173 

0.001916 

± 

0.000127 

-19 

0.004961 

± 

0.000197 

0.004387 

± 

0.000185 

0.002110 

± 

0.000467 

-18.5 

0.006454 

± 

0.000225 

0.007382 

± 

0.000838 

0.004008 

± 

0.000644 

-18 

0.007849 

± 

0.000869 

0.008353 

± 

0.000892 

0.005485 

± 

0.000753 

-17.5 

0.010007 

± 

0.000981 

0.012432 

± 

0.001088 

0.007384 

± 

0.000874 

-17 

0.012560 

± 

0.001099 

0.012238 

± 

0.001079 

0.016030 

± 

0.001095 

-16.5 

0.001432 

± 

0.001173 

0.012821 

± 

0.001105 

0.010337 

± 

0.001034 

-16 

0.017660 

± 

0.001303 

0.015599 

± 

0.001216 

0.013510 

± 

0.000850 

-15.5 

0.018052 

± 

0.001317 

0.014625 

± 

0.000971 

- 



-15 

0.030077 

± 

0.001505 

- 



- 



-14.5 

0.033572 

± 

0.001251 

— 



— 




z ~ 3 are noticeably higher than derived in nearly all pre¬ 
vious studies, while our inferred values of M* are therefore 
(unsurprisingly) somewhat lower. Our results in fact agree 
best with those recently derived by Weisz et al. (2014), who 
attempted to reconstruct the form of the UV LF out to 
z ~ 5 from the properties (including star-formation histo¬ 
ries) of the galaxies in the local group. We speculate that 
this agreement perhaps reflects the fact that our own study 
and that undertaken by Weisz et al. (2014) are the only 
studies to date which have probed to the depths required 
to properly determine the faint-end slope (in fact Weisz et 
al. (2014) reach down to Muv — —5 at z ~ 0), with an in¬ 


evitable resulting impact on the inferred best-fitting values 
of the other two Schechter parameters. While this agreement 
is interesting, and arguably impressive, we note that, unsur¬ 
prisingly, the uncertainties in our parameter values are much 
smaller than those presented by Weisz et al. (2014). 

Interestingly, the evolution of (fi* and M* derived here 
over the redshift range z ~ 2 — 4, is also very similar to that 
recently derived for the evolving emission-line galaxy LFs by 
Khostovan et al. (2015) (although their results were derived 
by locking the value of the faint-end slope, due to the lack 
of data of sufficient depth to constrain it). 
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Figure 8. Our derived best-fitting Schechter function parameter 
values at z ~ 2, 3, and 4, with the associated l-tr and 2-cr single¬ 
parameter confidence regions (corresponding to Ax^ = 1)4 after 
minimizing over the other parameter). As can be seen from the 
upper panel, a remains shallow, steepening little if at all over the 
redshift range. The lower panel shows that 0* drops gradually 
from z ~ 2 to z ~ 3, but then falls by a factor ~ 2.5 by z ~ 
4. However, this drop in number density is offset by a steady 
brightening in M* by ~ 1 mag. from z ~ 2 to z ~ 4, to the extent 
that UV luminosity density peaks at z ~ 2.5 — 3; see Fig. 11. 

5.2 Evolution up to z ~ 8 

In Fig. 10 we also attempt to place our findings in the wider 
context of the results derived from a number of recent stud¬ 
ies of the UV LF extending out to z ~ 8. It can be seen 
that the level of agreement is in fact better at z ~ 4 than at 
z ~ 2 or z ~ 3. The solid black lines in each panel of Fig. 10 
show simple parametric fits to the data (i.e. to the published 


Figure 9. The rest-frame UV (1500A) galaxy luminosity func¬ 
tions presented in Fig. 7, now overlaid to show the form of the 
evolution from z ~ 2 to z ~ 4. As is also clear from the Schechter 
function parameter values plotted in Fig. 8, it can be seen that 
i) the faint-end slope is little changed over this redshift range, ii) 
normalization drops only slightly between z ~ 2 and z ~ 3 but 
then drops by a factor ~ 2.5 by z ~ 4, and iii) the LF brightens 
steadily by ~ 1 mag. from z ~ 2 to z ~ 4. 

Schechter parameter values), to illustrate the overall evolu¬ 
tionary trend in each parameter as would be derived from 
the literature. These curves are meant to guide the eye, and 
are not meant to indicate our best estimate of true parame¬ 
ter evolution; indeed our new, more accurate determinations 
of 0* and M*atz~2 —3 clearly differ significantly from the 
literature average (for the well-understood reasons described 
above). Nevertheless, the evolutionary trend in a indicated 
by the simple straight-line fit shown in the bottom panel 
serves to re-emphasise how well our derived shallow faint- 
end slopes agree with the Weisz et al. (2014) results, and 
also shows that such values are in fact in very reasonable 
agreement with the general trend of a gradual steepening 
from a ~ —1.2 at z ~ 0 to a ~ —2 by z ~ 8. 


5.3 Luminosity Density 

We finish by considering the evolution of UV luminosity den¬ 
sity, puv, inferred from our LF determinations over the key 
redshift regime z ~ 2 — 4. While a full determination of star- 
formation rate density evolution also requires accounting for 
the substantial impact of dust obscuration, the luminosity- 
weighted integral of the UV LF does provide an important 
measurement of the unobscured star-formation rate density 
at each redshift. 

The results of this calculation are presented in Table 4 
and plotted in Fig. 11. Here we have deliberately performed 
all calculations using the LFs determined at Xreat = 1500 A 
to enable unbiased comparison of the derived values at each 
redshift. Table 4 includes results calculated at z ~ 2, 3 and 
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Hathi et al. 2010 
Sawicki & Thompson 2006 
van der Burg et al. 2010 
Bouwens et al. 2014 
Bouwens et al. 2007 
Bouwens et al. 2011 
Yoshida et al. 2006 
FInkelsteIn et al. 2014 
McLure et al. 2009,2013 
Ishigaki et al. 2014 
Atek et al. 2014 
Iwata et al. 2007 
Schenker et al. 2013 
Ouchi et al. 2009 
Schmidt et al. 2014 
Bradley et al. 2012 
Oesch et al. 2007 
Arnouts et al. 2005 
Wyder et al. 2005 
Reddy & Steidel 2009 
Alavi et al. 2014 
Sawicki 2012 
Weisz et al. 2014 
Oesch et al. 2010 
Bowler et al. 2014 
This Work 


z 

Figure 10. A compilation of the derived Schechter function parameter values for the UV galaxy LF over the redshift range z ~ 0 to 
z ~ 8, placing the new results derived in this paper, and the other results tabulated in Table 3, into the wider context of virtually all 
of cosmic time. Our own results, with error bars (see Fig. 8), are shown by the large black points, with other results from the literature 
plotted as indicated in the legend. The solid black line in each panel is a simple parametric fit to the data, plotted to illustrate the 
overall evolutionary trend in each parameter as inferred from the literature. The evolution of (p* and M* from z ~ 2 to z ~ 4 deduced in 
the present study is somewhat more dramatic than seen in most pre-existing direct studies of the LF at these redshifts, in fact agreeing 
best with the values inferred by Weisz et al. (2014) from galactic archaeology of the local group. The bottom panel again re-emphasises 
how well our derived shallow faint-end slopes agree with the Weisz et al. (2014) results, but also shows that such values are in fact in 
reasonable agreement with the general trend of a gradual steepening from a ~ —1.2 at z ~ 0 to a ~ —2 by z ~ 8. 


© 2014 RAS, MNRAS 000, [^?? 


















The galaxy UV luminosity function at z ~ 2 — 4 13 


Table 3. A compilation of the derived Schechter function parameter values for the UV galaxy LF over the redshift range z ~ 1.7 to 
2 ~ 4.0, placing the new results derived in this paper in the context of results presented in the literature over the last 10 years. The 
values tabulated here are included in Fig. 10 


Source 

z 

^rest / 

M* 

r 

a 

This Work 

1.7 

1500 

-19.61 ± 0.07 

0.00681 ± 0.00081 

-1.33 ±0.03 

Oesch et al. (2010) 

1.7 

1500 

-20.17 ±0.34 

0.00234 ± 0.00096 

-1.60 ±0.21 

Sawicki & Thomson (2006) 

1.7 

1700 

-19.80 ±0.32 

0.01698 ± 0.00489 

-0.81 ± 0.21 

Hath! et al. (2010) 

1.7 

1500 

-19.43 ±0.36 

0.00217 ± 0.00077 

-1.27 ± 0.00 

This Work 

1.9 

1500 

-19.68 ±0.05 

0.00702 ± 0.00066 

-1.32 ±0.03 

Oesch et al. (2010) 

1.9 

1500 

-20.16 ±0.52 

0.00219 ±0.00123 

-1.60 ±0.51 

Arnouts et al. (2005) 

2.0 

1500 

-20.33 ±0.50 

0.00265 ± 0.00020 

-1.49 ±0.24 

Alavi et al. (2014) 

2.0 

1500 

-20.01 ± 0.24 

0.00288 ± 0.00084 

-1.74 ± 0.08 

Weisz et al. (2014) 

2.0 

1500 

-19.36 ±0.28 

0.00750 ± 0.00400 

-1.36 ±0.11 

Hathi et al. (2010) 

2.1 

1500 

-20.39 ±0.64 

0.00157 ± 0.00115 

-1.17 ± 0.40 

Sawicki & Thomson (2006) 

2.2 

1700 

-20.60 ±0.38 

0.00301 ± 0.00176 

-1.20 ±0.24 

Sawicki (2012) 

2.2 

1700 

-21.00 ±0.50 

0.00274 ± 0.00024 

-1.47 ± 0.24 

This Work 

2.25 

1500 

-19.71 ± 0.07 

0.00759 ± 0.00088 

-1.26 ±0.04 

This Work 

2.25 

1700 

-19.99 ±0.08 

0.00620 ± 0.00077 

-1.31 ± 0.04 

Reddy & Steidel (2009) 

2.3 

1700 

-20.70 ±0.11 

0.00275 ± 0.00054 

-1.73 ±0.07 

Arnouts et al. (2005) 

2.7 

1500 

-21.08 ±0.45 

0.00162 ± 0.00090 

-1.47 ± 0.21 

Hathi et al. (2010) 

2.7 

1500 

-20.94 ±0.53 

0.00154 ± 0.00114 

-1.52 ±0.29 

This Work 

2.8 

1500 

-20.20 ±0.07 

0.00532 ± 0.00060 

-1.31 ± 0.04 

Arnouts et al. (2005) 

3.0 

1500 

-21.07 ±0.15 

0.00140 ± 0.00000 

-1.60 ±0.13 

Weisz et al. (2014) 

3.0 

1500 

-20.45 ± 0.26 

0.00410 ± 0.00200 

-1.36 ±0.13 

Sawicki & Thomson (2006) 

3.0 

1700 

-20.90 ±0.22 

0.00167 ± 0.00013 

-1.43 ±0.17 

van der Burg et al. (2010) 

3.0 

1600 

-20.94 ±0.14 

0.00179 ± 0.00051 

-1.65 ±0.12 

Reddy & Steidel (2009) 

3.05 

1700 

-20.97 ± 0.14 

0.00171 ± 0.00053 

-1.73 ±0.13 

This Work 

3.8 

1500 

-20.71 ± 0.10 

0.00206 ± 0.00033 

-1.43 ±0.04 

Bouwens et al. (2014) 

3.8 

1600 

-20.88 ±0.08 

0.00197 ± 0.00034 

-1.64 ± 0.04 

Bouwens et al. (2007) 

3.8 

1600 

-20.98 ±0.10 

0.00130 ± 0.00020 

-1.73 ±0.05 

Weisz et al. (2014) 

4.0 

1500 

-20.89 ±0.11 

0.00182 ±0.00010 

-1.58 ±0.08 

Sawicki & Thomson (2006) 

4.0 

1700 

-21.00 ±0.40 

0.00085 ± 0.00021 

-1.26 ±0.40 

van der Burg et al. (2010) 

4.0 

1600 

-20.84 ± 0.09 

0.00136 ± 0.00023 

-1.56 ±0.08 

Finkelstein et al. (2015) 

4.0 

1500 

-20.73 ±0.09 

0.00141 ± 0.00021 

-1.56 ±0.06 

Yoshida et al. (2006) 

4.0 

1500 

-21.14 ±0.14 

0.00146 ± 0.00041 

-1.82 ±0.09 


4, and also provides results for narrower redshift bins at 
2 ~ 1.7 (1.5 < 2 < 2.0) and 2 ~ 2.25 (2.0 < 2 < 2.5) to 
enable checking of the trend within the 2 ~ 2 redshift bin. 
We also provide results integrated down to different limiting 
absolute magnitudes: M 1500 = —17.7, M 1500 = —15, and 
A 7 i 5 oo = —10. The results to M 1500 = —17.7 are given for 
ease of comparison with many existing studies, while the 
convergence seen at the deeper limits shows that, because 
our derived faint-end slopes are fairly fiat, relatively little 
additional luminosity density is contributed by the faintest 
galaxies; Fig. 11 shows that puv has essentially converged 
by Misoo — —15. Regardless of the chosen integration limit, 
our results indicate that UV luminosity density (and hence 
unobscured star-formation rate density) peaks at 2 ~ 2.5—3, 
when the Universe was ~ 2.5 Gyr old. 

The formal uncertainties indicated by the error bars 
in Fig. 11 are fairly small, both because the integral of the 
LF is better constrained than the (somewhat degenerate) 
Schechter parameters, and because, with such a flat faint- 
end slope, uncertainties in a only have a minor effect on the 
luminosity-weighted integral. In practice, therefore, the true 
uncertainties are likely to be dominated by cosmic variance 
(although luminosity density is clearly less affected by cos¬ 
mic variance uncertainties than, for example, bright galaxy 
number counts). 


Finally, in Fig. 12 we show our derived UV luminosity 
densities (integrated down to M 1500 = —17.7) in the con¬ 
text of other recent determinations at comparable redshifts, 
and recent measurements extending to 2 ~ 9. Our new re¬ 
sults are more accurate than previous determinations, but 
in generally good agreement with existing results at 2 ~ 3 
and 2 ~ 4. At 2 ~ 2 our new result lies at the low end of 
the (widely discrepant) previously reported measurements, 
but is in fact still higher than the recent estimate provided 
by Alavi et al. (2014). Even allowing for cosmic variance, 
the basic conclusion that UV luminosity density peaks at 
2 ~ 2.5 — 3 appears secure. 

Of course, the precise epoch at which cosmic star- 
formation rate density reached a peak depends on the evo¬ 
lution of the correction for dust obscuration. At 2 ~ 2 a 
number of arguments indicate that this correction involves 
scaling the raw UV luminosity density by a factor of ~ 4 — 5 
(e.g. Reddy & Steidel 2009; Burgarella et al. 2013; Madau 
& Dickinson 2014), but whether this correction evolves sig¬ 
nificantly between 2 ~ 3 and 2 ~ 2 remains as yet unclear. 

Recent reviews of cosmic star-formation history based 
on data compilations have generally favoured a peak in cos¬ 
mic star-formation rate density at 2 ~ 2 (e.g. Behroozi, 
Weschler & Conroy 2013; Madau & Dickinson 2014) but at 
least some recent studies (e.g. radio: Karim et al. 2011; far- 
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Table 4. The rest-frame UV (1500A) luminosity densities as derived from our UV LFs from 2 ~ 1.7 to 2 ~ 4, with the luminosity- 
weighted integral performed down to three different magnitude limits. Because our derived faint-end slopes are fairly flat, relatively 
little additional luminosity density is contributed by the faintest galaxies, and Fig. 11 shows that puy has essentially converged by 
-^^500 — —15. Regardless of the chosen integration limit, it seems clear that UV luminosity density (and hence unobscured star-formation 
density) peaks at 2 ~ 2.5 — 3, when the Universe was ~ 2.5 Gyr old. The values given here are plotted in Fig. 11. 


2 

puv/^0^^ ergs s ^ Hz ^ Mpc ^ 

— 10 

pn^;/10^® ergss ^Hz ^Mpc ^ 

— 15 

puu/10^® ergss ^Hz ^Mpc ^ 

— 17.7 

1.7 

2.79 ±0.04 

2.62 ± 0.05 

1.91 ± 0.04 

1.9 

3.01 ± 0.04 

2.84 ± 0.06 

2.11 ± 0.07 

2.25 

3.13 ±0.04 

3.00 ± 0.06 

2.33 ±0.07 

2.8 

3.64 ± 0.01 

3.50 ± 0.04 

2.86 ±0.07 

3.8 

2.70 ±0.03 

2.56 ±0.05 

2.09 ±0.07 





Reddy&Steidel 2009,M=-17.5 
Bouwens et al. 2007 
Alavi et al, 2014, M=-17.7 
SawIckl&Thomson 2006 
McLeod et al, 2015 
This Work M=-17.7 


2 


2 


Figure 11. The rest-frame UV (1500A) luminosity densities 
as derived from our UV LFs from 2 ~ 1.7 to 2 ~ 4, with 
the luminosity-weighted integral performed down to three dif¬ 
ferent absolute magnitude limits: M1500 = —17.7 (green circles), 
-V/1500 = —15 (blue stars), and M1500 = —10 (pink triangles). 
The results to M1500 = —17.7 are shown for ease of comparison 
with many existing studies, while the convergence seen at the 
deeper limits shows that, because our derived faint-end slopes 
are relatively flat, relatively little additional luminosity density 
is contributed by the faintest galaxies (i.e. the plot shows that 
Puv has essentially converged by M 1500 — —15). Regardless of 
the chosen integration limit, it seems clear that UV luminosity 
density (and hence unobscured star-formation density) peaks at 
2 ~ 2.5 — 3, when the Universe was ~ 2.5 Gyr old. The values 
plotted here are tabulated in Table 4. 


infrared: Burgarella et al. 2013; emission-line: Khostovan et 
al. 2015) favour a peak nearer 2 3. The latter scenario 

is more obviously consistent with the new UV results pre¬ 
sented here, but a definitive answer awaits more direct mea¬ 
surements of dust-enshrouded star-formation at 2 2 — 4 

from forthcoming deep sub-mm/mm surveys with SCUBA-2 
on the JCMT and the Atacama Large Millimeter Array. 


Figure 12. Our derived UV (1500A) luminosity density val¬ 
ues at 2 ~ 2, 3 and 4, with integration performed down to 
Mi 5 oo = —17.7, are here plotted as the black points, and com¬ 
pared to results of similar calculations performed by other au¬ 
thors as indicated in the legend. Our results are more accurate 
than previous determinations, but in generally good agreement 
with existing results at 2 2 :: 3 and 2 ~ 4. At 2 ~ 2 our new result 
lies at the low end of the (widely discrepant) previously reported 
measurements. 


6 CONCLUSION 

We have exploited the high dynamic range provided 
by combining the Hubble Ultra Deep Field (HUDF), 
CANDELS/GOODS-South, and UltraVISTA/COSMOS 
surveys to derive a new, robust measurement of the evolv¬ 
ing rest-frame ultraviolet galaxy luminosity function over 
the key redshift range from 2 2 to 2 4. 

The unparalleled multi-frequency photometry available 
in this survey ‘wedding cake’, combined with the (relative) 
wealth of deep optical and near-infrared spectroscopy in 
these fields, has enabled us to derive accurate photometric 
redshifts for c::: 95% of the galaxies in the combined survey 
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(with a reliability and accuracy that are competitive with 
the very best achieved to date, as verified in Appendix A). 

This has then enabled us to assemble robust and com¬ 
plete galaxy samples within redshift slices at z ~ 2, 3 and 4, 
facilitating a new determination of the form and evolution of 
the UV galaxy LF, that probes ~ 3 — 4 magnitudes fainter 
than previous (unlensed) surveys at z ~ 2 — 3, and does 
not rely on potentially incomplete colour-colour selection 
techniques. The SED fitting undertaken to determine the 
photometric redshifts has also allowed us to determine ac¬ 
curate rest-frame UV absolute magnitudes (Misoo or Miroo, 
as required for comparison with previous results). 

Our new determinations of the UV LF extend from 
Misoo — —22 (AB mag) down to Misoo = —14.5, —15.5 and 
— 16 at z ~ 2, 3 and 4 respectively. Fitting a Schechter func¬ 
tion to the LF data as determined from the Vmax estima¬ 
tor, at 2 ~ 2 — 3 reveals a much shallower faint-end slope 
(a = —1.32 ± 0.03) than the steeper values (a ~ —1.7) re¬ 
ported by Reddy & Steidel (2009) or by Alavi et al. (2014) 
(who utilised gravitional lensing to help sample the faint 
end of the LF). By performing the Schechter function fit¬ 
ting down to differing limiting magnitudes, we show that our 
measurement of faint-end slope is robust (i.e. the inferred 
value plateaus/converges at Misoo > —17). By 2 ~ 4 the 
faint-end slope has steepened slightly, to a = —1.43 ± 0.04. 
Although these values are significantly shallower than the 
aforementioned pre-existing estimates, we find they are in 
excellent agreement with the values recently inferred by 
Weisz et al. (2014) (from galactic archaeology of the local 
group), and are in fact consistent with the overall evolution¬ 
ary trend in a from 2 = 0 to 2 = 8, as gleaned from a review 
of the literature. 

Analysis of the other best-fitting Schechter function pa¬ 
rameters reveals that our derived number density normaliza¬ 
tion, (j >*, is higher than nearly all previous estimates at 2 ~ 2 
(except, again, Weisz et al. 2014), declines only slightly by 
2 ~ 3, and then drops by a factor ~ 2.5 to 2 ~ 4 (where 
our value agrees well with most previous measurements). 
Meanwhile, this drop in number density is offset by a steady 
brightening in M* by ~ 1 mag. from 2 ~ 2 to 2 ~ 4, to the 
extent that UV luminosity density does not drop signifi¬ 
cantly until the negative density evolution takes over and 
dominates beyond 2 ~ 3. 

Finally, we have compared our new UV LF determina¬ 
tions, and the resulting inferred evolution of UV luminosity 
density {puv), with results from a range of previous studies 
extending from 2 ~ 0 out to 2 ~ 9. Because our new mea¬ 
surements yield fairly flat faint-end slopes, our estimates of 
puv are relatively robust; they have essentially converged 
by Muv — —15, and are little influenced by remaining un¬ 
certainties in a. We conclude that unobscured UV luminos¬ 
ity density (and hence unobscured star-formation density) 
peaks at 2 ~ 2.5 — 3, when the Universe was ~ 2.5 Gyr old. 
Whether or not this coincides with the peak in total cosmic 
star-formation rate density (psfr) depends on the results of 
ongoing efforts to determine the level and evolution of dust 
obscuration at this epochs. 
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APPENDIX A: PHOTOMETRIC REDSHIFT 
ACCURACY AND RELIABILITY 


In this appendix we provide some additional details on the 
reliability and accuracy of our photometric redshifts in the 
three survey fields utilised in this study, and assess how our 
results compare with other recently published photometric 
redshift catalogues. 

To estimate the accuracy of the photometric-redshift 
estimation procedure, we compare our photometric redshifts 
with their spectroscopic counterparts, for the subsamples of 
galaxies for which high-quality spectroscopic redshifts are 
known. 

Following standard practice, we use the following statis¬ 
tics to quantify the accuracy and reliability of the photomet¬ 
ric redshifts. 

First, the basic scatter, cr, around the Zphot '■ Zspec line 
is defined as: 

a = rms[Az /(1 + Zspec)] (Al) 

where Az = Zphot — Zspec- 

Second, ‘catastrophic outliers’ are defined as galaxies 
for which: 


\Az\/{l + Zsp^c) > 0.15. 


(A2) 


Third, the scatter can be recalculated after exclusion of 
the catastrophic outliers (in order to estimate the tightness 
of the core Zphot ■ Zspsc relation); this measure of scatter is 
usually denoted as as- 

Fourth, an alternative measure of scatter, that mini¬ 
mizes the impact of (but does not require the removal of) 
catastrophic outliers, is the normalised median absolute de¬ 
viation of Az defined as: 


aNMAD = 


1.48 X median{ 


|A^| 

1 Zspec 


(A3) 


As described in Section 3, we have assembled sub¬ 
samples of galaxies with high-quality spectroscopic redshifts 
within each of the three fields. In the HUDF field there 
are 218 such galaxies, and we obtained acceptable pho¬ 
tometric redshifts (i.e. < 50) for 210 of these. The 

Zphot : Zspsc plot for this sub-sample of 210 galaxies is shown 
in the upper panel of Fig. Al, with the normalized redshift 
error, Az/{1 -\- Zspec), plotted against Zspec shown in the 
lower panel. The values of the four aforementioned statis¬ 
tics are given in the upper panel. Analogous plots are then 
shown for the corresponding subsamples of 2677 galaxies in 
CANDELS/GOODS-S (Fig.A2) and 1671 galaxies in Ul- 
traVISTA/COSMOS (Fig. A3). In all three fields the spec¬ 
troscopic redhsifts provide coverage from z ~ 0 to at least 
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Figure Al. Zpf^^t versus z^pec for the 210 galaxies in the HUDF 
field with high-quality spectroscopic redshifts and acceptable pho¬ 
tometric redshifts < 50). The outlier/scatter measurements 
are given in the upper panel, while the lower panel redshift error 
as a function of z^p^c, with the dashed blue lines inidcating the 
0.15 boundary used to define the catastrophic outliers. 

z ~ 5, and the outlier/scatter statistics are consistent and 
competitive with the accuracy of the very best photomet¬ 
ric redshifts as reported elsewhere in the recent literature; 
outlier fraction is always significantly lower than 5%, and 
i^NMAD — 0.026 in all three fields. 

Within the CANDELS/GOODS-S field, an alternative 
set of photometric redshifts has recently been released by the 
3D-HST team (Skelton et al. 2014). In Fig. A4 we plot our 
own Zphot '■ Zspec results for this field (for the same 2677 
galaxies shown in Fig. A2) along with the corresponding 
results as derived from the 3D-HST photometric redshift 
catalogue, and in Table Al we compare the resulting out¬ 
lier/scatter statistics. Clearly these two photometric redshift 
catalogues are of comparably high quality, although we note 
that the outlier fraction achieved here is significantly lower, 
possibly because the 3D-HST photometric catalogue does 


Figure A2. ZpjjQ^ versus z^pec for the 2677 galaxies in the 
CANDELS/GOODS-S field with high-quality spectroscopic red¬ 
shifts and acceptable photometric redshifts (y^ < 50). The out¬ 
lier/scatter measurements are given in the upper panel, while the 
lower panel redshift error as a function of z^pec, with the dashed 
blue lines inidcating the 0.15 boundary used to define the catas¬ 
trophic outliers. 

not contain the HST E-band imaging. The slightly smaller 
i^NMAD achieved by the 3D-HST team appears to result 
from improved accuracy at low redshifts, possibly driven by 
their inclusion of medium-band ground-based Subaru imag¬ 
ing. However, at the redshifts of interest in the present study 
(z > 1.5), our own measurements yield a slightly smaller 
scatter than is achieved by using the 3D-HST catalogue. 

Similarly, within the HUDF field, an alternative set of 
photometric redshifts has recently been released by Rafelski 
et al. (2015). In Fig. A5 we plot our own Zphot '■ Zspec results 
for this field along with the corresponding results as derived 
from the Rafelski et al. (2015) photometric redshift cata¬ 
logue. Here we are plotting results for 207 galaxies (because 
3 of the 210 galaxies plotted in Fig. Al do not have photo¬ 
metric redshifts in the Rafelski et al. (2015) catalogue). We 
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Figure A3. 2p^ot versus Zspec for th® 1671 galaxies in the Ul- 
traVISTA/COSMOS field with high-quality spectroscopic red- 
shifts and acceptable photometric redshifts (x^ < 20). The out¬ 
lier/scatter measurements are given in the upper panel, while the 
lower panel redshift error as a function of Zspea with the dashed 
blue lines inidcating the 0.15 boundary used to define the catas¬ 
trophic outliers. 


also include the results from both of the alternative redshift 
estimation techniques used by Rafelski et al. (2015), which 
are based on the Bayesian Photometric Redshift (BPZ; Ben¬ 
itez 2000) and Easy and Accurate Zphot from Yale (EAZY; 
Brammer, van Dokkum &; Coppi 2008) algorithms. The for¬ 
mer uses a set of PEGASE SED models which have been 
re-calibrated based on observed photometry and spectro¬ 
scopic redshifts from the FIREWORKS catalogue; emission 
lines are included and luminosity functions observed in COS¬ 
MOS, GOODS-MUSIC and the UDF are used as priors. The 
latter method uses the default EAZY SEDs, with emission 
lines again included. In Table A2 we again compare the re¬ 
sulting outlier/scatter statistics. Despite the fact that the 
Rafelski et al. (2015) results include new WFC3/UVIS pho¬ 
tometry, again it can be seen that all three sets of results 



Figure A4. A comparison of the reliability/accuracy of our pho¬ 
tometric redshifts in the CANDELS/GOODS-S field with that 
achieved by the 3D-HST team (Skelton et al. 2014). Our results 
(red points) are overlaid on the results derived from the public 
catalogue released by the 3D-HST team (green points) for the 
same spectroscopic sub-sample of 2677 galaxies as previously dis¬ 
cussed and presented in Section 3 and Fig.A2. Outlier fractions 
and scatter statistics are summarized in Table Al, and discussed 
in the text. 


Table Al. The reliability and accuracy of the photometric red¬ 
shifts for galaxies in the CANDELS/GOODS-S field, as achieved 
here (see Section 3) and alternatively by utilising the the public 
3D-HST photometric redshift catalogue (Skelton et al. 2014). cr, 
^NMAD ^rid font (see text) have been calculated for the high- 
quality spectroscopic sub-sample of 2677 galaxies (see Fig. A4). 
Our own results yield a lower outlier fraction, possibly because 
the 3D-HST photometric catalogue does not contain the HST Y- 
band imaging. The slightly smaller achieved by the 3D- 

HST team appears to result from improved accuracy at low red¬ 
shifts, possibly driven by their inclusion of medium-band ground- 
based Subaru imaging. However, at the redshifts of interest in the 
present study {z > 1.5), our own measurements yield a slightly 
smaller scatter than is achieved by using the 3D-HST catalogue. 
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font 

This Work 

0.120 

0.036 

0.027 

3.0% 

3D-HST 

0.135 

0.027 

0.013 

4.3% 


are competitive (presumably because our own calculations 
utilise the VIMOS f/-band ground-based photometry, mini¬ 
mizing the additional impact of the new UVIS data). Indeed, 
our outlier fraction is lower than yielded by the Rafelski et 
al. (2015) EAZY results which, as in our own calculations, 
avoid the use of luminosity function priors, and our cata¬ 
logue yields the lowest value of cjnmaD’ 
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The galaxy UV luminosity function at z ~ 2 



Figure A5. A comparison of the reliability/accuracy of our pho¬ 
tometric redshifts in the HUDF field with that achieved by the 
Rafelski et al. (2015). Our results (red points) are overlaid on the 
results derived using the BPZ (green triangles) and EAZY (blue 
triangles) algorithms by Rafelski et al. (2015) for 207 of the 210 
galaxies in the HUDF spectroscopic sub-sample as previously dis¬ 
cussed and presented in Section 3 and plotted in Fig. Al. Outlier 
fractions and scatter statistics are summarized in Table A2, and 
discussed in the text. 


Table A2. The reliability and accuracy of the photometric red- 
shifts for galaxies in the HUDF, as achieved here (see Section 
3) and alternatively by utilising the the new public photometric 
redshifts released by Rafelski et al. (2015). a, crjsiMAD and font 
(see text) have been calculated for 207 of the galaxies with high- 
quality spectroscopic redshift in the HUDF (see Fig. A5). Our 
own results produce a lower outlier fraction than yielded by the 
Rafelski et al. (2015) EAZY results which, as in our own calcula¬ 
tions, avoid the use of luminosity function priors. Moreover, our 
catalogue yields the lowest value of cTTVAfAD- 
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